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Abstract 

The study of protostellar collapse or infall has proven 
unusually difficult and controversial. Some historical con¬ 
text will be provided, against which recent advances can 
be measured. We now have a substantial number of ob¬ 
jects with signatures that can be interpreted in terms of 
collapse, but a number of issues remain. One issue is the 
effect of chemical variations, especially depletion in the 
dense, cold interiors of cores that are likely to form low 
mass stars. Strategies for dealing with this issue depend 
on using dust emission to constrain the density and tem¬ 
perature distribution, leaving molecular line observations 
to constrain dynamics and abundance distributions. Re¬ 
cent progress in this area will be described and we will 
consider the next challenges to be overcome. Interferome- 
teric observations, especially with future instruments, will 
provide a powerful tool. In combination with chemical 
studies coupled with dynamical models, the observations 
made possible by interferometers should hnally put this 
subject on a firm foundation. 
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1. Introduction 

While gravitational collapse is the accepted explanation 
for star formation, the observational evidence for it is 
very scarce. Establishing this picture observationally is 
a strong reason for studying infall. This subject has long 
been controversial (cf. Goldreich & Kwan 1974 and Zuck- 
erman & Evans 1974), and the burden of proof has rested 
heavily on those who interpret observations in terms of 
collapse. While it will always be difficult to rule out other 
explanations for observational signatures of collapse, the 
collapse interpretation is beginning to beneht from the 
principle of simplicity: other explanations are becoming 
increasingly contrived. 


If we can establish the basic picture, we can move on 
to testing particular theories for how the collapse pro¬ 
ceeds. There are many theoretical models that differ in 
their predictions of the density and velocity fields. In prin¬ 
ciple, these predictions may be subject to observational 
tests. As an example, collapse according to the Larson- 
Penston similarity solution (Larson 1969; Pension 1969) 
produces much wider linewidths than inside-out collapse 
from a singular isothermal sphere (Shu 1977). The line 
profiles expected for the Larson-Penston solution do not 
match data from regions forming low mass stars (Zhou 
1992) but may find application in regions forming mas¬ 
sive stars. 

Finally, observational constraints on the dynamics will 
provide an alternative way to assess the timescales for 
star formation. At present, these are determined only by 
statistical arguments based on the number of objects in 
different observational classes, assumed to be an evolu¬ 
tionary sequence. This method is subject to a number of 
uncertainties, including selection biases and the possible 
importance of factors other than age, such as environ¬ 
ment, etc. For example, the usual view is that protostars 
spend 2 x 10^ yr in the Class 0 phase and 2 x 10® yr in 
the Class I phase. The former number is based on the 
fact that the p Ophiuchus cloud has ten times as many 
Class I objects as Class 0 objects (Andre & Montmerle 
1994). However, Visser, Richer, & Chandler (2001, 2002) 
have found as many Class 0 objects as Class I objects 
in a recent survey of other regions, implying a timescale 
for the Class 0 stage that is 10 times longer than usu¬ 
ally assumed. Direct observation of infall speeds in var¬ 
ious objects would provide an independent estimate of 
timescales. 

The topic of infall has been reviewed by Zhou & Evans 
(1994), Myers (1997), Zhou (1999), Evans (1999), and 
Myers, Evans, & Ohashi (2000). In this paper, I intend 
mostly to update those reviews, but I will begin with 
some background. 


Proceedings of the conference “Chemistry as a Diagnostic of Star Formation,” University of Waterloo, Canada, 21-23 
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2. Dark Ages and Renaissance 


Partly because of the historical controversy referred to 
above, the study of collapse suffered through the “dark 
ages” from the mid-1970s to the mid-1990s. During this 
period, every claim to find evidence for collapse was fol¬ 
lowed by counter-claims with alternative explanations [cf. 
Goldreich & Kwan (1974), Liszt et al. (1974) vs. Zuck- 
erman & Evans (1974); Snell & Loren (1977) vs. Leung 
& Brown (1977); Walker et al. (1986) vs. Mundy, Myers, 
& Wilking (1986) and Menten et al. (1987)]. For further 
historical discussion and references, see Evans (1991) and 
Myers et al. (2000). The net effect of all this contention 
was to make the study of infall seem disreputable. 

In addition to sociological factors, physical factors 
have made the field difficult. Collapse motions tend to 
be small compared to other motions and hence easily ob¬ 
scured. A rough estimate of the infall speed is given by 
the following: 


Vinf ~ 


1 km s ^ 


- M,/Mq i 0-5 
-r/lOOOAU. 


( 1 ) 


Thus, for low-mass stars, one needs to probe small 
radii to see substantial infall motion. Other motions are 
often substantially larger. For a sample of low-mass cores, 
Caselli & Myers (1995) found that the non-thermal (tur¬ 
bulent) contribution to the linewidth (FWHM) could be 
fitted by the following relation: 

Autur-b = 1-3 km (2) 


While this relation does not necessarily apply to the 
linewidth as a function of radius for a given region (see 
Goodman et al. 1998), there is some evidence that tur¬ 
bulent linewidths are less at small radii on average. Thus 
improved resolution is an obvious route to improve the de¬ 
tectability of infall. However, the decreases in linewidth 
will be limited as thermal broadening takes over. The 
minimum possible linewidth (aside from narrowing by 
maser emission) is 


Autb = 


8(ln2)fcTif 


= 0.215 km s 


-1 


Tk 


m(amu) ’ 


( 3 ) 


where Tk is the gas kinetic temperature and m is the 
mass of the species whose line is being observed. Based 
on this equation, one would gain by observing heavier 
molecules, but purely thermal broadening is not usually 
seen in heavier molecules. In fact, turbulence often per¬ 
sists at the subsonic level even on small scales. The sonic 
linewidth is 

^Vsomc = = 0.14 km (4) 

y '^mean 


where rrimean is the mean mass per particle. In some very 
quiescent clouds, the turbulence is considerably less than 
this value. 

In addition, rotational motions will increase as regions 
closer to the forming star are probed. If angular momen¬ 
tum is conserved, 

Wrot tx 7’“^ (5) 

Ohashi et al. (1997) found that angular momentum ap¬ 
pears to be conserved inside about 6000 AU. On much 
smaller scales, where the circumstellar disk forms, angu¬ 
lar momentum will be redistributed, but the rotational 
velocity does appear to increase with decreasing radius 
in the infall region. 

Finally, and worst, outflows seem to begin very soon 
after a central source has formed, and they are very en¬ 
ergetic at early times. These move matter at velocities 
ranging up to a few times 100 km s“^ in the wind and 
a few times 10 km s“^ in the molecular outflow. Most 
troubling for the study of infall is the low velocity mo¬ 
tion of dense gas very near the base of the outflow. These 
motions can easily obscure or confuse infall signatures. 

The renaissance in the study of infall can be traced to 
four developments. First, the discovery of objects in very 
early stages provided more fertile grounds for infall stud¬ 
ies. Presaged by extensive surveys for dense gas in nearby 
dark clouds by Myers and colleagues (e.g., Myers, Linke 
& Benson 1983), the identification of Class 0 objects by 
Andre et al. (1993) and Class —1 objects (also called Pre- 
Protostellar Cores or PPCs) by Ward-Thompson et al. 
(1994) was a crucial development. Second, simple models 
for the density and velocity fields during collapse were 
presented by Shu (1977), but the idea of inside-out col¬ 
lapse from an initial, singular isothermal sphere really 
came to prominence after the seminal review article by 
Shu, Adams, feLizano (1987). Third, Zhou (1992) con¬ 
verted the theoretical models of both the inside-out col¬ 
lapse and the Larson-Penston solution into simulated line 
profiles as a function of time, clarifying the characteris¬ 
tic line shapes expected for collapse. Finally, a credible 
example of collapse was found in the small globule B335 
(Zhou et al. 1993). 

The key signature of infall is a line profile that ideally 
is double-peaked with a central self-absorption dip, with 
the blue peak being stronger than the red peak (Leung 
& Brown 1977, Zhou 1992, Myers et al. 1996). Cartoons 
that demonstrate why this signature arises can be found 
in Figs. 1^ and This signature will appear only if the 
line has a suitable optical depth and critical density. Ob¬ 
servations of an optically thin line that peaks near the ve¬ 
locity of the self-absorption dip are necessary to rule out 











Evans: Studying Infall 


3 


two separate velocity components. Some optically thick 
lines will not show a self-absorption dip, but will show a 
blue-skewed profile. The latter are referred to by Myers 
et al. (1996) as “red-shoulder” profiles, but they are the 
same as the “blue-skewed” profiles. For simplicity, I will 
refer to either of these as “blue profiles.” 



Antenna 



Figure 1. The origin of various parts of the line profile for 
a cloud undergoing inside-out collapse. The static enve¬ 
lope outside Vinf produces the central self-absorption dip, 
the blue peak comes from the back of the cloud, and the 
red peak from the front of the cloud. The faster collapse 
near the center produces line wings, but these are usually 
confused by outflow wings. 


The full double-peaked signature was seen clearly in 
several lines toward B335 by Zhou et al. (1993) and mod¬ 
eled in more detail by Choi et al. (1995). While some 
issues have arisen later regarding B335 (§4), it was credi¬ 
ble enough to return some respectability to the field, and 
it led to the next development. 

3. Surveys for Infall Signatures 

One example is insufficient to verify a theory for low- 
mass star formation, so a number of surveys began for 
blue profiles. These initially used the molecular tracers 
that showed blue profiles in B335, CS and H 2 CO. Various 
lines were used for the optically thin tracer. In a survey 
of 12 globules, similar to B335, Wang et al. (1995) found 



Figure 2. A schematic explanation of why line profiles 
of optically thick, high-excitation lines are skewed to the 
blue in a collapsing cloud. The ovals are loci of con¬ 
stant line-of-sight velocity, for v{r) oc Each line of 

sight intersects these loci at two points. The point closer 
to the center will have a higher Tgx, especially in lines 
that are hard to excite, so that Tex{R 2 ) > Tex(i?i) and 
Tex{B 2 ) > T'ex(Bi). If the line is sufficiently opague, the 
point i?i will obscure the brighter R 2 , but B 2 lies in front 
of Bi. The result is a profile with the blue peak stronger 
than the red peak (Zhou & Evans 1994). 


3 blue profiles using H 2 CO as the optically thick tracer. 
Mardones et al. (1997) surveyed a large number of Class 
0 and I sources, finding 14 out of 37 blue profiles in CS 
and 15 out of 47 in H 2 CO. Most surprisingly, Lee et al. 
(1999) found blue profiles in 17/70 starless cores, using 
CS. In the inside-out picture, a central luminous object 
forms almost as soon as collapse begins, so a phase with 
collapse and no central luminosity would be very rare. 

Meanwhile, it turned out that HCO'*' showed even 
deeper self reversals (Fig. ||) than CS and H 2 CO toward 
B335, and searches commenced using HCO'*' J = 3 ^ 2 
as the optically thick line and H^^CO'^ as the optically 
thin line. Gregersen et al. (1997) found blue profiles to¬ 
ward 6 out of 18 Class 0 sources. Because the study by 
Mardones et al. (1997) using CS and H 2 CO had found 
fewer blue profiles among the Class I sources than among 















4 


Evans: Studying Infall 


Class 0 sources, it was important to put the Class O/I 
searches on a common footing. A study of Class I sources 
by Gregersen et al. (2000) found blue profiles in 8 out of 
16 Class I sources and Gregersen & Evans (2000, 2001) 
found blue profiles in 6 of 17 Class —1 cores (starless cores 
with submillimeter emission). Clearly the phenomenon 
of blue profiles was widespread, but far from universal. 
Worse yet, some lines showed red profiles! 



4 6 8 10 12 

V (km/s) 

Figure 3. The line profile of the J = 3 ^ 2 HCO^ line 
toward B335 (Evans, Rawlings, & Choi 2002b). 

To quantify the results, Mardones et al. (1997) intro¬ 
duced the excess: E = {Nuue - Nred)/Ntot, where Nuue 
and Nred are the number of blue and red profiles in the 
total sample of A^tot objects. A line was counted as blue 
if the peak of the optically thick line was shifted blue- 
ward of the peak of the optically thin line by a quarter 
of the linewidth of the optically thin line: 6v = {vthick — 
Vthin)/^Vthin < —0.25. Similarly, a red line would have 


Excess versus 

Class. 



Class 

-1 

0 

I 

Excess 

0.30 

0.31 

0.31 


Sv > 0.25. The excess should be zero on average if the 
skewing of line profiles were caused by random outflow 
blobs or other effects without some overall inward mo¬ 
tion. 

The surveys did reveal positive excesses, suggesting a 
systematic tendency toward inward motions. One ques¬ 
tion that arose was whether there was any tendency for 
the excess to vary with time. In some pictures, the main 
infall occurs during the Class 0 stage, with much less in¬ 
fall occurring during the Class I stage (e.g., Andre et al. 
2000). Studies with CS and H 2 CO did find a lower value 
of E among Class I sources (Mardones et al. 1997), but 
studies with HCO+ found the same E = 0.31 in Class I 
sources as in Class 0 sources (Gregersen et al. 2000). Be¬ 
cause HCO+ is more opaque than CS or H 2 CO, it may 
be that it was better able to reveal infall at later stages 
when less material remains in the envelope. Using HCO'*", 
a surprisingly large excess was found for Class —1 cores 
by Gregersen & Evans (2000), but the sample was very 
small. The value in Table 1 is based on counting only Sv 
that are less than —0.25 by at least one cr; if all those 
with Sv < —0.25 are counted, E = 0.50. 

Another approach taken with surveys is to ask if dif¬ 
ferent tracers agree. Gregersen et al. (2000) assembled a 
table of 40 sources, most with data in CS, H 2 CO, HCO''", 
and HCN. There were some cases in which the tracers 
disagreed, but considerably more in which they agreed. 
In quite a few cases, the HCO+ showed blue profiles, 
while the others showed no asymmetry. Gregersen et al. 
(2000) suggested that the HCO+ remained optically thick 
enough to develop a blue profile at later stages, when the 
GS and H 2 CO were no longer opaque enough. They pre¬ 
sented some simple simulations of line profiles to support 
this suggestion. A better understanding of chemical ef¬ 
fects in these regions is needed to put this suggestion on 
firmer footing. 

4. Complications 

While the surveys indicated a substantial excess of blue 
profiles over red profiles, the presence of red profiles raises 
a red flag. Without understanding those profiles, how can 
we be sure what the blue profiles mean? The red profiles 
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are likely to be caused by dense clumps in outflowing gas. 
Clearly some of the blue profiles could be caused by the 
same thing. While the statistics suggest that the predom¬ 
inance of blue prohles indicate collapse, a blue profile in 
an individual core may not be a reliable signature of col¬ 
lapse. 

In view of the troubled history, modern students of in¬ 
fall have generally been quite cautious, subjecting a blue 
profile to a number of tests before describing it as a cred¬ 
ible example of collapse. This process reminds me of the 
classic work of Puritan allegory (Bunyan 1688), so I offer 
a modern allegory for the “Profile’s Progress” in Fig. A 
blue profile could just represent the presence of two veloc¬ 
ity components, which must be ruled out by observations 
of an optically thin line that peaks near the dip in the op¬ 
tically thick line, rather than showing two peaks. If there 
is an embedded source (Class 0 or I), a map of the source 
should reveal a peak on or near the source; the peak could 
be in the “blueness,” measured in various ways, or in the 
depth of the dip. The peak of the blue component need 
not be on the source if rotation is also present, but a pat¬ 
tern consistent with expectations for rotation should be 
seen. Now, the hardest test appears; a blue profile must 
resist the temptation to just be an outflow blob! The pro¬ 
file that passes this test may be considered a candidate 
for collapse. Agreement with someone’s detailed model 
is required to achieve the status of a credible example. 
Opinions differ as to whether any regions have achieved 
this status. 

If no embedded source is known (Class —1), the path 
is slightly different (no outflow to worry about) and is 
not shown in Figure First, without a source, one has 
to check instead that some indicator is at least centered 
on the peak of the submillimeter continuum emission. In 
these cases, one may be seeing core formation, rather than 
star formation, and the pattern of blue profiles may be 
quite extended (e.g., Tafalla et al. 1998). The blue pro¬ 
files cannot be tracing infall onto a central object, so they 
must be compared to models of core formation by am- 
bipolar diffusion or turbulent dissipation (e.g., Myers & 
Lazarian 1998) or to models of collapse before a central 
object forms (Lee et al. 2003). 

Because infall velocities should increase at smaller 
radii, it seemed that interferometers might provide clearer 
evidence of infall. Indeed, some studies did look promis¬ 
ing (Chandler & Sargent 1993, Zhou et al. 1996, Wilner 
et al. 1997). However, Choi et al. (1999) found some dis¬ 
crepancies. In particular, the red peak was very weak or 
even absent in some cases. This was an extreme exam¬ 
ple of a trend already seen in single-dish data toward 


A Profile’s Progress 



Figure 4- The progress of a blue profile through the many 
pitfalls on the path toward “salvation, ” as a credible ex¬ 
ample of collapse (with apologies to John Bunyan). 

higher blue/red peak ratios in the data than predicted 
by the models (Gregersen et al. 1997). In addition, in¬ 
terferometers can artificially produce what appear to be 
self-absorption dips by resolved out the extended emis¬ 
sion from a static or slowly moving envelope (e.g., Gueth 
& Guilloteau 1999). 

When the CS J = 5 ^ 4 line was imaged in B335 by 
Wilner et al. (2000), they found that the line prohle, while 
blue, was narrower than predicted from the best htting 
model to the single dish data (Choi et al. 1995). At a 
minimum, this result suggested lower velocities at small 
radii than predicted from the inside-out collapse model. 
Even worse, the map showed that the blue emission was 
displaced from the infrared source and thus more likely 
an outflow blob. 

Around the same time, it was becoming clear that 
many molecules can be heavily depleted in the dense, cold 
regions of Class —1 sources (e.g., Caselli et al. 1999), and 
that these depletions might affect the line profiles strongly 
(e.g., Rawlings and Yates 2001) if they persist into the 
infalling stage. 

The detection of blue profiles in many Class — 1 cores 
also called into question the basic model of inside-out 
collapse, suggesting instead that inward motions of some 
kind are commonly seen before formation of a central 
luminosity source. Maps of the lines indicated extended 
regions of blue prohles, inconsistent with very early stages 
of an inside-out collapse. (Tafalla et al. 1998; Lee, Myers, 
& Tafalla 1999, 2001) 
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All these developments suggested the need for a more 
robust tracer of the density structure than provided by 
molecular lines because detailed simulation of the line 
profiles required models of the density, temperature, ve¬ 
locity, and abundance structure. 

5. Using Dust and Gas 

While dust emission or absorption provides no direct in¬ 
formation on velocities, studies of dust can be quite useful 
for the infall problem. Studies of dust continuum emission 
or absorption can constrain the distribution of density 
and temperature, cutting down the number of free param¬ 
eters in models of line profiles. Conveniently, several new 
instruments became available in the late 1990s that rev¬ 
olutionized observations of dust emission. For low-mass, 
cold regions, SCUBA on the JCMT, operating at 850 
and 450 /xm, and various bolometer arrays on the IRAM 
30-m dish have provided many maps of submillimeter 
or millimeter emission (e.g., Shirley et al. 2000; Ward- 
Thompson et al. 1999; Visser et al. 2001, 2002). 

Studies of cores through absorption (or scattering) by 
dust are complementary to emission studies. These have 
a long heritage, but one has to use absorption in the in¬ 
frared to probe the denser regions. Studies in the near- 
infrared by Lada et al. (1994) and in the mid-infrared 
using ISO (Bacmann et al. 2000) are examples of this 
technique. The use of the mid-infrared will be extended 
by observations with SIRTF of dense cores (Evans et al. 
2002a). 

Proper interpretation of the emission data requires 
knowledge of the dust temperature distribution (Tdir)) 
because at the low temperatures around low luminos¬ 
ity sources, the dust emission is quite sensitive to tem¬ 
perature. Unlike the case of luminous sources, the emis¬ 
sion is not in the Rayleigh-Jeans limit and the Td(r) 
is not a power law (Shirley et al. 2000); consequently, 
convenient analytical approaches are inappropriate. For¬ 
tunately, there are publicly available codes to compute 
Td(r) self-consistently for a given density distribution, in¬ 
cluding both internal and external heating (Egan, Leung, 
& Spagna 1988, Nenkova, Ivezic, & Elitzur 2000). 

For line profile modeling, one needs to know TK(r), 
the distribution of gas kinetic temperature. For much of 
the relevant range of radii, the density is high enough that 
collisions of molecules with dust grains makes Tj^ ~ T^, 
but in the lower density envelope, a separate energetics 
calculation is needed to obtain the kinetic temperature. 
While one might worry that heavy depletion of molecular 
coolants deep in the interiors might affect Tk, Goldsmith 


(2001) has shown that the effects are small. Once the den¬ 
sity is high enough to affect the coolants significantly, it is 
high enough to enforce Tk ~ Td to reasonable accuracy. 

The results of models with self-consistent calculations 
of Td{r) (Leung 1975, Evans et al. 2001) show that Class 
— 1 cores are very cold in the center {Td ~ 7 K) and 
warmer on the outside {Td ~ 14 K if exposed to a typical 
interstellar radiation field). Density distributions that are 
constant in the center but approach power laws at large 
radii, such as Bonnor-Ebert spheres fit the radial intensity 
profile well (Ward-Thompson, et al. 1999, Evans et al. 
2001). The low central Td decreases the emission from the 
center, even at millimeter and submillimeter wavelengths. 
When the lower temperatures in the center are accounted 
for, the central densities of Class —1 cores with reasonably 
strong submillimeter emission range from 3 x 10^ to 10® 
cm“® (Evans et al. 2001). 

At the very low temperatures and high densities deep 
in Class —1 cores, substantial depletion of the molecules 
onto dust grains is expected on relatively short timescales. 
The observational evidence clearly confirms this expecta¬ 
tion (e.g., Rawlings 2000, Caselli 2002 for reviews). Many 
Class —1 cores are nearly invisible in the usual molecular 
tracers of column density, such as C^®0. Other species 
that form from less condensible precursors or that profit 
from the depletion of their primary destroyers follow the 
dust emission more closely (e.g., Lee et al. 2003). Espe¬ 
cially favored are molecules derived from N 2 , ions, and 
deuterated species. On these grounds, N 2 D+ may be the 
last gaseous molecular tracer still working after the others 
have chilled out (Caselli et al. 2002). 

This depletion makes it very hard to trace the ve¬ 
locity field in the interior, but studies at larger radii 
show widespread regions of blue profiles, as described ear¬ 
lier. These profiles can be successfully modeled in some 
cases with a Plummer-like model (Whitworth & Ward- 
Thompson 2001); in particular, Lee et al. (2003) have fit 
the HCO+ line profiles toward L1544 at several positions 
with such a model. The density and temperature distribu¬ 
tions of this model are nearly identical to the best fitting 
Bonnor-Ebert sphere and thus predict dust emission that 
matches observations, but the Plummer-like model also 
provides a velocity distribution needed to compare to line 
profiles. In the case that fits L1544, the inward velocity 
peaks around 0.015 pc. An alternative approach (Bourke 
et al. 2002) to the same core treats each molecular tran¬ 
sition with a two-layer model; when applied to different 
molecular transitions, one infers a velocity field larger at 
smaller radii, qualitatively consistent with the modeling 
by Lee et al. (2003). The exciting thing about these re- 
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suits is that we are beginning to probe velocity fields in 
very early stages. These will set timescales and initial 
conditions for the formation of central objects. However, 
these are only initial steps and much work remains, even 
to show that these motions are necessarily related to col¬ 
lapse. 

Once a central object forms, the core is heated from 
the inside as well as the outside, and Td rises in the core. 
The submillimeter emission from dust becomes stronger 
and easier to map. Analysis of a substantial sample of 
Class 0 and I objects (Shirley et al. 2002; Young et al. 
2003) finds that the radial intensity profile of the submil¬ 
limeter emission is fitted well by power laws in density 
[n(r) = nf{r/rf)~P]. An interesting dichotomy appears: 
for sources that are nearly spherical, {p) = 1.8; for the 
small number of sources that are quite elongated, typical 
p ~ 1. While the mean value of p is similar to that of 
the singular isothermal sphere, these objects should have 
evolved into the inside-out collapse stage, which predicts 
p < 1.5 inside ri„/. No sign of a shift to a smaller value of 
p at small radii is apparent in the data from single dishes, 
which probe r > 2000 AU, in Class 0 sources (Shirley et 
al. 2002). More surprisingly, most Class I sources also 
show no sign of a change in p outside 2000 AU (Young 
et al. 2003), even though Vinf should be larger in these 
sources. 

For B335, the inside-out collapse model that fits the 
molecular line data (Choi et al. 1995) does not fit the dust 
continuum data, which prefer a power law density. Studies 
of extinction in the near-infrared using HST observations 
(Harvey et al. 2001) find a power law with p = 2 in the 
outer parts of B335, with some evidence for a turnover 
near the found by Choi et al. (1995). The extinction 
method runs out of stars inside r ~ 3500 AU, but recent 
studies of dust emission with the IRAM interferometer 
probing scales from 500 to 5000 AU were fitted well by a 
power law withp = 1.65±0.05 (Harvey et al. 2002). These 
results suggest that the inside-out model does reproduce 
the density distribution, returning the issue to the models 
of molecular lines. 

Combining the constraints from dust continuum emis¬ 
sion with those of molecular lines in a single, coherent 
model is clearly necessary. Hogerheijde & Sandell (2000) 
demonstrated the power of this technique; they found 
good fits to models of inside-out collapse for three of four 
sources. The fourth source, L1489 IRS, was later modeled 
as a more evolved object with a rotating, contracting disk 
(Hogerheijde 2001). Such studies are clearly necessary to 
reduce the large number of free parameters once molecu¬ 
lar depletion is considered. 


6. The Role of Interferometers 

While data from interferometers must be interpreted with 
care (see discussion in §4), they provide the only probe 
of the velocity fields on small scales. Choi (2002) has ex¬ 
plored some of the systematic effects and simulated line 
profiles from infall as seen by interferometers such as the 
SMA. As interferometers acquire more antennae, their 
imaging fidelity will improve, but it will remain impor¬ 
tant to account for extended emission. 

In addition to looking for infall signatures in optically 
thick lines, interferometers have been used to probe ve¬ 
locity fields in optically thinner tracers (e.g., Ohashi et al. 
1996, Ohashi 1997, Momose et al. 1998). The signature 
of infall is a velocity gradient along the projected minor 
axis of a flattened structure that is larger than the disk. 
Rotation produces a velocity gradient along the projected 
major axis, allowing one to separate these two motions. 
Typically, both infall and rotation are inferred from the 
data. As always, confusion with the outflow is a problem, 
but this technique is likely to prove powerful as we begin 
to understand which species do not trace the outflow. 

In addition to studies of emission, interferometers pro¬ 
vide a new possibility. Because they are sensitive to emis¬ 
sion from disks, the continuum emission from the disk can 
serve as a background lamp and, with even higher spatial 
resolution and sensitivity, an occulting disk for the back 
half of the cloud (Evans 2001). For disk emission that is 
sufficiently bright, the front of the cloud will appear in 
absorption while the rest of the cloud appears in emission, 
producing an inverse P-Cygni profile. With instruments 
like ALMA, one may be able to resolve the opaque part of 
the disk, occulting the back of the cloud altogether, leav¬ 
ing only the absorption line from the front of the cloud. A 
redshift in this absorption line is the “smoking gun” for 
collapse. Hints of this absorption were seen toward NGC 
1333 IRAS 4A by Choi et al. (1999), but the IRAM ob¬ 
servations of Di Francesco et al. (2001) show a much more 
dramatic inverse P-Cygni profile. The remaining question 
in this source is whether the absorption arises from the 
infalling envelope or from a foreground cloud that just 
happens to be redshifted. Whatever the outcome in this 
particular case, this technique will be routine with in¬ 
struments like ALMA and will allow detailed study of 
the velocity field (Evans 2001). 

7. Studies of the Velocity Eield 

For the moment, most studies of the velocity field are still 
done with single dishes. By matching line profiles to mod¬ 
els over a map of a source or by modeling the line profiles 
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toward the center for lines with different excitation and 
opacity, it is possible to decode the velocity field, albeit 
with ample uncertainties, in part due to unknown deple¬ 
tion profiles. There are two general approaches to this 
problem. One is to take a well-defined theoretical model 
and try to constrain it with many observations to deter¬ 
mine whether it fits the data with any reasonable set of 
input parameters (e.g., Zhou et al. 1993, Choi et al. 1995, 
Lee et al. 2003). The other approach is to decode the ve¬ 
locity field empirically by treating each molecular line in 
a fairly simple way, but using the ensemble to constrain 
v^nfir)■ 

The latter approach forms the basis for a study of 
IRAM 04191 by Belloche et al. (2002); this paper repre¬ 
sents the current state of the art. This source appears 
to be in a very early stage of infall. The line profiles 
of optically thick, hard-to-excite lines toward the center 
show blue profiles, with the self-absorption dip moving 
farther redward as the excitation requirements of the line 
increase. By excitation requirements, I mean some com¬ 
bination of temperature and density. Belloche et al. inter¬ 
pret these line profiles in terms of a constant Vmf =0.1 
km s“^ for r > 3000 AU, and Vinf increasing inward for 
r < 3000 AU. From line profiles away from the center, 
they find evidence of rotation as well, with Vrot increasing 
with decreasing radius in to r = 3500 AU, but constant 
(with large uncertainties) inside 3500 AU. This source 
is destined to become one of the standard test cases for 
infall studies, along with B335 and L1544. 


8. Future Prospects 

While many questions remain and skepticism is justified, 
the field seems to be moving from arguing about infall to 
studying it. The combined constraints of dust and molec¬ 
ular line data will be important, as will a deeper the¬ 
oretical understanding of the depletion of molecules and 
other chemical effects. Future study of velocity fields, just 
getting under way, will require self-consistent modeling 
tools, especially 2D or 3D radiative transfer. Observa¬ 
tions will need to probe a wide range of spatial scales, 
and theory will have to provide testable models that can 
be generated by observers from publicly available codes 
to compare against their data. If we acquire a better un¬ 
derstanding of infall in regions forming low mass stars, 
we may be in a position to address the question in the 
much more complex regions forming massive stars. 
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